Shell corrections for finite depth potentials with bound states only 
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A new method of calculating unique values of ground-state shell corrections for finite depth 
potentials is shown, which makes use of bound states only. It is based on (i) a general formulation 
of extracting the smooth part from any fluctuating quantity proposed by Strutinsky and Ivanjuk, 
(ii) a generalized Strutinsky smoothing condition suggested recently by Vertse et al., and (iii) the 
technique of the Lanczos a factors. Numerical results for some spherical heavy nuclei ( 132,154 Sn, 
i80,208pk anc j 298 are presented and compared to those obtained with the Green's function 
oscillator expansion method. 
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The binding energies of nuclei as a function of nu- 
cleon numbers and shape parameters show an irregular 
behaviour superimposed on a smooth trend. This irreg- 
ularity (shell corrections) is mainly due to the gross-shell 
non-uniformity of the single-particle (sp) levels near the 
Fermi energy. The smooth trend can be reproduced well 
by the liquid-drop or droplet mass formula Q, while the 
shell corrections can be calculated with the Strutinsky 
method [2j.y(. 

The Strutinsky method converges very well (there is an 
interval of values 4] of both the energy smoothing param- 
eter and the order of the curvature-correction polynomial 
where the shell correction practically remains the same) 
when it is applied to infinite depth potentials like the har- 
monic oscillator potential or a Nilsson model potential. 
Difficulties appeared 0|, however, in calculations using 
finite depth potentials where the number of bound levels 
above the Fermi energy was not enough for the appli- 
cation of the traditional energy-averaging (iJ-averaging) 
Strutinsky procedure 0, Q • 

In several papers, e.g., Refs. [1,|EI3> ideas were dis- 
cussed to improve the convergence of the ^-averaging 
procedure for finite potentials by including the contin- 
uum effect on the level density. Other methods H, 0, 
have also been applied to obtain shell corrections of finite 
potentials. Very recently in Ref. the authors have 

drawn the conclusion that the proper treatment of the 
continuum effect on level density does not guarantee that 
the standard -E-averaging procedure converges. Conse- 
quently, an alternative prescription for defining shell cor- 
rections of finite potentials, based on the assumption of a 
local linear energy dependence of the smoothed level den- 
sity (generalized Strutinsky smoothing condition), was 
suggested. Moreover, a practical method of calculating 
the smooth sp level density in combination with the pro- 
posed definition of shell correction was used in Ref. 



to obtain unambiguous shell corrections of the finite de- 
formed Woods-Saxon potential for drip-line nuclei. This 
new procedure was called Green's Function Oscillator Ex- 
pansion (GFOE) method. 

Apart from the methods mentioned above, there are 
prescriptions based on avera ging in the nucleon num- 
ber space (TV-averaging) [II H El E E3- The A- 
averaging method does not need any unbound levels. The 
difference between the A-averaging and the ^-averaging 
procedures is due to the symmetry of the sp hamiltonian, 
i.e., the degeneracy of the sp levels ^3- The A- and E- 
averaging are particular cases of a general formulation of 
extracting the smooth part from any fluctuatin g qu antity, 
suggested by Strutinsky and Ivanjuk in Ref 13]. Un- 
fortunately, the versions of the A-averaging pre scrip tion 
discussed so far in the literature (e.g., Refs. [l3Lll4l llq) 
do not seem to yield stable results. 

In this paper a new method of calculating the ground- 
state shell corrections of finite depth potentials is pro- 
posed, which makes use of bound states only. The 
method is based on (i) the general formulation shown 
in Ref. |13| . (ii) the generalized Strutinsky smoothing 
condition suggested recently in Refs. El[l2]> and (iii) 
the technique of the Lanczos a factors jlq . 

Method: (i) In the general formulation presented in Ref. 
|13| (when it is referred to the nucleon number space, A- 
averaging) the smooth part E s (A) of the total sp energy 
sum E S (N) (A is the nucleon number of the system ) 
was determined as the least-square-deviation (LSD) fit 
to EJN) by an M-th degree polynomial of A 1 / 3 [HQ 
Q,El- The smoothed total sp energy E S (N) is expressed 
as 



E S (N) = Cm(N,X)E s (X), 



(1) 
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where Ai and A2 are the lower and the upper limits of 
the averaging interval Ai < A < A2. From now on we 
will fix Ai and A2 to the minimum and the maximum 



2 



nucleon number that the potential- well can contain, i.e., 
Ni = 1 and N 2 = N max which is the sum of the de- 
generacy of all individual bound states. It is important 
to point out that the degeneracy of a sp level will not 
be taken into account in Q with a degeneracy factor, 
because the nucleons numbered by A, belonging to a de- 
generated sp level, appear with different weights w(X) in 
the smoothing function £m which reads as 



Cm(N, X) = K m (N, X)w(X). 



(2) 



The curvature-correction Km{N,X) is the polynomial 
of the M-th degree, which is composed of polynomials 
Pk (X) orthonormal in the discrete interval 1 < A < N max 
with the weight w(X) [l9| and reads as 



K M (N, X) =J2Pk(N) Pk (X), 



(3) 



fe=0 



where Y^\=i° Pk(X)p k ' (A)w(A) = S kk ' . The smoothed 
total sp energy E S (N) in expression Q can be written 
in terms of the sp energies e\ and the smoothed sp occu- 
pation numbers n\ as 



E S (N) 



£ 

A=l 



(4) 



41 A- 1 / 3 MeV being the mean distance between the gross- 
shells. Consequently, we follow the generalized Strutin- 
sky smoothing condition suggested in Refs. 0, 0] to 
obtain an optimal value of the parameter A (A op ) for 
a given M. This condition was supported by extensive 
calculations 0, of the smoothed sp level density us- 
ing both the semiclassical Wigner-Kirkwood method and 
the exact expression of the level density which incorpo- 
rates the continuum effects properly. From these studies, 
it was observed that the smoothed level density shows 
a linear energy dependence in the intermediate energy 
region of bound states. In Refs. the following 

generalized smoothing condition was suggested: In an 
energy interval [e; , e„] which is wider than the mean dis- 
tance huiQ between the gross-shells (e.g., e u — t\ = 2hwo), 
the deviation of the smoothed sp level density g(e, A, M) 
from linearity should be minimal. Hence, A op for a given 
M is calculated by minimizing 



X 2 (A,M) 



[g(e, A,M) - a- befde, (6) 



where the parameters a and b are uniquely determined 
for each value of A and M by the method of least squares. 
The smoothed sp level density g(e, A, M) is obtained as 
follows: The smoothed total sp energy E S (X) (A refers 
to the number of particles with which the sp levels are 
filled, i.e. 1 < A < N max ) can be written in terms of g(e) 
as 



where n x = J2v=\* (m(N,v) and Y,\=F 71 a = N. The 
conservation of the number of particles is guaranteed by 
the definition of the smoothing function ( M The 
argument of the smoothing function Qm is understood 
as x\ rather than A, where the variable x\ is a linear 
function of A 1 / 3 , i.e., x x = (A 1 / 3 - iV 1/3 )/A. The param- 
eters Nq and A determine the maximum and the width 
of the weight function w(x\) which will be a Gaussian 
w(x\) = ex.p(—x\)/y/Tr centered at the Fermi level, i.e., 
the parameter Nq is fixed at the nucleon number TV for 
which the shell correction in the iV-averaging procedure 



SU(N) = E S (N) - E S (N) 



(5) 



E S (X) 



eg(e)de, 



(7) 



with fi\ determined by the condition of conservation of 
the number of particles 



){e)de. 



(8) 



From these two equations, g(fi\) and fi\ can be calcu- 
lated as 



/Oa) = ( 



d 2 E s (X), l 



dX 2 



(9) 



is calculated. Therefore, the remaining free parameters 
of the equations above are A and the M degree of the 
smoothing function An unambiguous (converged) 

value of the shell correction 8U(N) means that it should 
not depend (or very weakly at most) on A and M, which 
will be discussed below. 

(ii) In our calculations we have not found any plateau 
of 5U(N) regarding the parameter A, for a fixed M 
degree of the smoothing function (m- The value of 
A is searched around A = 1, which approximately 
corresponds to the energy smoothing parameter (7 = 
l.2hujo) of the traditional .E-averaging procedure, ftwo = 



dE s (X) 
dX : 



(10) 



where E S (X) are obtained with the moving average (QJ. 
The expressions (l§))- (|10f) allow the construction of the 
function g = g(e, A,M). The energy interval [e;,e u ] is 
centered around the half of the energy of the lowest sp 
level (i.e., ~ 0.5ei). The search for A op begins at a small 
A value below A = 1, e.g., A = 0.5 where the effect of the 
shell fluctuations on g(e, A, M) is still present. A is then 
gradually increased until the first minimum in x 2 (^? 
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appears at A op . This A op represents the smallest value of 
A for a given M that smooths out the shell fluctuations 
effect on g(e, A, M). The shell correction 5U(N, A op , M) 
is the optimal one for a given M. This prescription is 
repeated for higher values of M. If variations of 8U with 
M are small (e.g., < 0.2 MeV), then the mean value of 
those SU represents a unique value of the shell correction. 
The improvement of the convergence with respect to M 
will be discussed next. 

(iii) Numerical instabilities happen concerning the con- 
vergence of the shell corrections with increasing M degree 
of the smoothing function £m. This is due to the effect 
of Gibbs oscillations shown by Lanczos in the theory of 
applied Fourier-analysis in Ref. [Taj , which always ap- 
pear when a smooth function is expanded in terms of a 
complete orthonormal set of functions and the expansion 
is truncated at a finite number N of terms. A technique 
to strongly decrease the instabilities and accelerate the 
convergence of the shell corrections may use the so-called 
Lanczos a factors 0, |2(], HJ- The main idea of this 
method consists of damping the higher order terms of 
the function expansion multiplying the expansion coef- 
ficients by attenuation factors (0 < < 1). The 
factor is 1 only for the first value of k, e.g., k — 0; 
for increasing k the decreases monotonously and be- 
come almost zero for the highest subscript k = N. The 
fc-dependence of is arbitrary and can be chosen to fit 
the actual problem best. 

The a factors enter the shell correction method dis- 
cussed above through expression @ for the curvature- 
correction polynomial Km- For a given nucleon number 
N, expression can be identified as a truncated expan- 
sion of the true curvature-correction function in terms of 
the orthonormal polynomials pk(X), where the expansion 
coefficients are pk(N). Therefore, the modified coeffi- 
cients &\Pk(N) will replace the coefficients pk(N). As 
in Refs. [2(| |2l| , the following of factors 

m l-exp{-[o(fc-M-l)/(Af+l)] 2 } 
1 — exp(— or) 

are used, where the value a = 5 was found in trials as in 
Refs. |2(j, |21| to optimize the convergence rate, i.e., the 
formation of a plateau. 

Numerical results: The bound states used in the 
present calculations are obtained with the Potential Sep- 
arable Expansion (PSE) method proposed by Revai in 
Ref. 201. The PSE's rigorous foundation was de- 
mostrated in Ref. pH ] and its efficiency was shown in 
several papers, e.g., Refs. |20j,|22|. For the sake of sim- 
plicity, the finite depth nuclear potential is chosen to be 
the following spherical Woods-Saxon with a spin-orbit 
term 

V(r) = -VoM + 1a(A)v I^ (1 . B)j (12) 

2 mc r ar 

where h/mc = 0.21 fm, / and f so are the same func- 
tion, but with different parameters, i.e., /( so )(r) = {1 + 



TABLE I: Potential parameters used to calculate the bound 
states: The potential depth is given in MeV, the radii and 
diffusenesses in fm, while the strength A of the spin-orbit po- 
tential is adimensional. See text for further details. 



Protons 



Nucleus 


Vb 


ro 


a 


A 


I SO 




Rc 






m gn 6 


59.94 


1.275 


0.70 


36 


1.30 


0.70 


6.62 


10 


7 


208p| 3 a 


59.21 


1.286 


0.762 


18.29 


0.88 


0.60 


6.69 


10 


7 


298 114 6 


59.62 


1.275 


0.70 


36 


1.30 


0.70 


8.86 


12 


10 


180p b i« 


53.39 


1.275 


0.70 


36 


1.30 


0.70 


7.34 


10 


7 








Neutrons 










132 Q b 

bn 


39.26 


1.347 


0.70 


36 


1.30 


0.70 




10 


7 


208p b c 


44 


1.27 


0.67 


32 


1.27 


0.67 




10 


7 


298 114 c 


43 


1.27 


0.67 


32 


1.27 


0.67 




12 


10 


154 Sn i, 


34.64 


1.347 


0.70 


36 


1.30 


0.70 




10 


7 
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exp[(r— ro( so ) J 4 1 / 3 )/a ( so )]} _1 . For protons, the Coulomb 
potential Vcoui is taken to be that of a uniformly charged 
sphere with charge (Z — l)e (Z being the total charge of 
the system) and the radius R c , which is added to expres- 
sion |Q . 

We will present numerical results for spherical heavy 
(stable) nuclei like 132 Sn, 208 Pb and 298 114 as well as for 
drip-line nuclei like 154 Sn and 180 Pb. The potential pa- 
rameters along with the maximal number of nodes n max 
and partial waves l m ax included in the harmonic oscil- 
lator basis for converged values of the sp energies are 
given in Table 1. In a realistic calculation, the poten- 
tial parameters should be selected in such a way that 
the experimental sp energies around the Fermi level are 
reproduced. 

Fig. 1 shows the optimal proton shell corrections 6U P 
for 208 Pb as a function of the M degree of the smoothing 
function Cm- The dependence of SU P on the a param- 
eter of the fife factors [a = 4 (circle), 5 (squares) and 6 
(triangles)] as well as the effect of these factors on the 
plateau formation of SU P (e.g., comparing squares with 
stars) can also be seen in Fig. 1. It is observed that the 
best plateau is formed for a = 5, and the Lanczos factors 
are hugely important to determine a unique value of the 
shell correction. As in Ref. it was found that the 

values of A op are strongly correlated with M, i.e., in the 
plateau of SU A op increases with increasing M. 

In Fig 2, we compare the smoothed sp level density 
of the neutrons in 208 Pb, which results from expressions 
Pl-iPfl for M = 10 at A£ p = 0.99 (solid curve), to 
that obtained by us with the GFOE method (dashed 
curve). The number of oscillator shells included in the 
basis for the GFOE-calculation is N osc = 36, the degree 
of the smoothing polynomial is 10, and the resulting op- 
timal energy smoothing parameter 7 op = 1.57 'tvuiQ. The 
smoothed sp level densities are very similar to each other 
in the middle of the energy region where their linearity 
is required by the generalized Strutinsky smoothing con- 
dition. The remaining oscillations of the solid curve in- 
dicate that the effect of the shell fluctuations has not yet 
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FIG. 1: Optimal proton shell correction 6U P for Pb as a 
function of the degree M of the smoothing function £m • The 
lines are to guide the eyes. See text for further details. 



TABLE II: Shell corrections (in MeV) obtained with the 
present approach 5U are compared to those SE calculated 
with the GFOE method 0]- See text for further details. 



Neutrons Protons 



Nucleus 


SU n 


SE n 


5U P 


SE P 


ly2 Sn 


-8.70 


-9.18 


-7.40 


-7.10 


208p b 


-11.74 


-11.00 


-8.85 


-8.12 


298 114 


-9.37 


-8.40 


-3.78 


-3.90 


154 Sn 


-1.20 


9.30 






180p b 






-0.22 


-8.0 



been completely eliminated by the found value of A™ . In 
contrast to the dashed curve, the solid curve abruptly in- 
creases near the continuum threshold and this behaviour 
is the same as that of the semiclassical level densit y ob - 
tained with the Wigner-Kirkwood method (see Ref . 11]). 
It is important to stress that unlike the i5-averaging pro- 
cedures such as the GFOE method or the semiclassical 
method, in the approach suggested in the present work, 
which is based on the iV-averaging, it is not the smoothed 
sp level density that determines the shell corrections, but 
only its linear behaviour in the intermediate energy re- 
gion. 



FIG. 2: The smoothed sp level density of the neutrons in 
208 Pb, resulting from the present approach (solid curve), is 
compared to that obtained by us with the GFOE method 
(dashed curve). See text for further details. 



Fig. 3 shows optimal neutron (upper part) and proton 
(lower part) shell corrections for 132 Sn (circles), 208 Pb 
(squares), 298 114 (triangles), 154 Sn (stars) and 180 Pb (di- 
amonds) as a function of the degree M of the smooth- 
ing function £m- A plateau is fairly well formed in all 
cases for M around 10. The heavier the system (i.e., 
more sp levels are included), the better the plateau can 
be seen. In table 2, the mean value of the optimal shell 
corrections belonging to the plateau SU are compared 
to the shell corrections SE obtained with the GFOE 
method. The shell corrections of the present approach 
can be larger (e.g., for neutrons in 132 Sn or for protons 
in 298 114) or smaller than those resulting from the GFOE 
method. The same features are observed comparing the 
shell corrections obtained with the semiclassical Wigner- 
Kirkwood method and the GFOE method (see Ref. lUl). 
For well-bound nuclei, the maximal deviation found in ta- 
ble 2 is for neutrons in 298 114 and is about 1 MeV. The 
difference cannot be explained by the so-called symmetry 
correction [TH Il6| which is always positive and was found 
in calculations to be about 2-3 MeV for the nuclei stud- 
ied. It seems that the generalized Strutinsky smoothing 
condition, used here in combination with the N-averaging 
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FIG. 3: Optimal neutron (upper part) SU„ and proton (lower 
part) SUp shell corrections for 132 Sn (circles), 208 Pb (squares), 
298 114 (triangles), 154 Sn (stars) and 180 Pb (diamonds) as a 
function of the degree M of the smoothing function (m- The 
lines are to guide the eyes. See text for further details. 



procedure, removes large part of the difference between 
the shell corrections arising from the iV-averaging and the 
-E-averaging prescriptions for well-bound nuclei. For 
the two drip-line nuclei studied, namely 154 Sn (neutrons) 
and 180 Pb (protons), it was found that the optimal shell 
corrections also shows a plateau with respect to M, but 
in this case the plateau largely deviates from the shell 
correction provided by the GFOE method (see table 2). 
Since reasons for such a discrepancy are still unknown, 
we do not recommend the use of the present method to 
drip-line nuclei for the time being. Further investigation 
is needed to clarify this point. 

Conclusions: To conclude we would say that ground- 
state shell corrections for finite depth potentials can be 
calculated using bound states only, by means of the 
new method proposed in the present paper. The re- 
aver/aging procedure to obtain the smoothed total sp en- 
ergy has been combined with both the generalized Struti- 
nsky smoothing condition and the Lanczos a factors. The 
new method provides unique values of the shell correc- 
tions and strongly reduces the well-known difference be- 
tween the iV-averaging and the i?-averaging prescriptions 
for well-bound (stable) nuclei, but this is not the case 
for drip-line nuclei. The method can be useful to cal- 
culate potential energy surfaces for fusion and fission in 
the macroscopic-microscopic approach with a two-center 
shell model based on two realistic Woods-Saxon poten- 
tials without the need of including unbound sp levels. 
Works in this direction are in progress. 
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